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I. INTRODUCTION 



During last years, in investigations of complex, self-organizing, fractal structures and 



various physical phenomena such as subdiffusion, turbulence and chemical reactions, as we 

social and biological systems, Tsallis [l 



as various economica 
Sharma-Mittal 



34 



19], Renyi 



20 



26 



26 



331, 



1 



35] statistics as well as superstatistics |36l441| are extensively used along 



with Gibbs one. The 



?sallis entropy is widely applied in various directions of nonextensive 



statistical mechanics 



^6 

n 



54. 



plasma 
tion of equi 
in Refs. 



13j. Some examples are the phenomena of subdiffusion 



42 



451 ] and turbulence |46r|53j, and the investigations of transport coefficients in gases and 



551 ] . as well as quantum dissipative systems 



56( . The problems of the construc- 



ibrium thermodynamics in the framework of generalized statistics were discussed 



ation 
Refs. 



68 



57H61]. The energy fluctuations 571]. kinetics of nonequilibrium plasma 62]. prob- 

the ozone layer 64j, universality in non-Debye relax- 
671 ] were investigated within the Tsallis formalism. In 



lems of self-gravitating systems 



and complex systems 



63 
66 



- 170] Tsallis statistics was applied to description of chemical reactions, in particu 



lar, nonlinear equations of react ion- diffusion processes were obtained in Ref. [681 ] . Despite 
the wide application of Tsallis entropy as a gene ralization of Gibbs-Shannon entropy, the 
Renyi entropy is of great interest as well 



71 



731 ] . In particular, in this case it is possible 



to determine a connection between t 
It is important to note the papers 174 



re index q and the heat capacity of the system [29] 



76 



78] by Luzzi, Vasconcellos and Ramos where the 



nonequilibrium statistical operator (NSO) method and the Renyi entropy are applied to 
description of systems far from equilibrium. In particuar, the nonequilibrium ^-dependent 
Renyi ensemble as well as the generalized distribution functions of bosons and fermions was 

physical-chemical systems 
75] . Therein investigations 



obtained [74J. A statistical approach for description of fracta 
based on non-Fickian diffusion processes was proposed in Ref. 
of anomalous diffusion in fractal-like electrodes in microbatteries were carried out. The 
experiments on anomalous luminescence at nanometer quantum dots in semiconductor het- 



erostructures were described 



76] in such approach. The Refs. 79 



8l] are devoted to investi- 



gation of nonlinear kinetics based on the Kramers, Boltzmann and Fokker-Planck equations 
within the framework of generalized statistics. 

In the present paper, one approach to formulation of an extensive statistical mechanics 



of nonequilibrium processes 33] is considered, based on the Zubarev NSO method 82 
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84| and the maximum entropy principle for the Renyi entropy. This statistical approach 
is applied to a consistent description of reaction-diffusion processes in the "gas-adsorbate- 
metal" system. React ion- diffusion and adsorption-desorption processes on the metal surface 
are nonlinear, they manifest oscillation character, possess memory effects and are actual in 
terms of nanostructures formation on surfaces in catalytic phenomena 85 



II. RENYI ENTROPY AND NONEQUILIBRIUM STATISTICAL OPERATOR 
METHOD 

The nonequilibrium state of a classical or quantum system of interacting particles is 
completely described by the nonequilibrium statistical operator g(x N ; t), which satisfies the 
classical or quantum Liouville equation: 

j t Q{x N -t)+iL Ne {x N -t) = 0. (1) 

iLjv is the Liouville operator of a system. 

Within the NSO method framework, we will be looking for solutions of Eq. (TjQ) which 
are independent of the initial conditions. The solutions will depend explicitly only on the 
observable quantities 

J dT N P nQ {x N -t) = {P n ) t . (2) 
The nonequilibrium statistical operator: 

g(x N ; t) = 6rel (x N ; t) - [ e £ ^T(t, t') (1 - P rd (t')) iL N g rel (x N ; t')dt', 



where T(t,t') = exp + |— f*, (1 — P re i(t')) iL^dt'j is the evolution operator containing the 
projection. The relevant statistical operator (distribution function) Qrei{ xN ]t) will be deter- 
mined using the maximum entropy principle for the Renyi entropy 

= in / dr N g%t) (3) 



1-q 

at fixed parameters of the reduced description, taking into account the normalization con 
dition. 

Relevant statistical operator corresponding to the Renyi entropy maximum: 



q-l 



drv 1 



g-l 



l 

9-1 



(5) 



Zji(t) is the partition function of the relevant statistical operator, 5P n = P n — (P n Y- The 
Lagrange multipliers F n (t) are defined from the self-consistency conditions: 



(PnY — (PnYreli 

where (. . = / dTV . . . Q re i(x N ; t). 

The generalized Kawasaki- Gunton projection operator: 

Sg re i(t) 



(6) 



Prel(t)g' = 



dr N P ne '. 



(7) 



An action of the operators P re i(t)iL N on the relevant statistical operator can be pre- 
sented as follows, by means of generalized projection which now acts on dynamic variables, 
P rel (t)iL N g rel (t) = P rel (t)A(t)g rel (t) = (P(t)A(t))g rel (t), 



p(t)... = (...) t rel + J £ s 



(...SP, 



m/rel ' 



(8) 



Here, f mn {t) is determined by the relations f mn (t) = ^Jf> = (|^y) , 



HPY 
5F(t) 



i - -{5Pr\t)Y rel F 



-{5Pr\t)5py rel = f{ti 



(9) 



(10) 



(11) 



where I is the unit matrix and 5 [...] = [...] — ([. . .])* ei . We also use the notation 

^ n 

Taking into account that (1 - P re i(t))iL N g re i(t) = -J2 n I n (t)F n (t)g re i(t), where 

J„(t) = (l-P(t))V 1 (t)> n 

are the generalized flows, we can now write down an explicit expression for the nonequilibrium 
statistical operator 

g(x N ;t) = g rel (x N ;t) + £ f e^'^T(t, t')I n {t')F n {t')g rel {x N ■ t')dt' . (12) 

This allows us to obtain the generalized transport equations for the reduced-description 
parameters. They can be presented in the form 

^{PrnY = (PmYrel + £ l[ WM')^')^, (13) 



with the generalized transport kernels (memory functions) 

tp mn (t,f) = J dT N {P m T(t,t')I n (t')g rel (t')} 
which describe the dissipative processes in the system. 



(14) 



III. REACTION-DIFFUSION PROCESSES 

Let us start with the Hamiltonian of the system of "gas-adsorbate-metal" in the following 
form H = H' + H reac , H' = H a + H l a nt . Here, H a is the Hamiltonian of the gas subsystem 
considered according to classical approach; H™ 1 is the Hamiltonian describing interaction 
between gas atoms and atoms adsorbed on the metal surface; H reac is the Hamiltonian 
of interaction for chemical reactions between atoms or molecules adsorbed on the metal 
surface 85] 



Hreac = ^ ((a', b'\$ reac \a, b)qf,q^qa,qb + («', b'\$ reac \a, b)*q£q£ q a 'qi') , (15) 

a,b,a',b' 

where (a', b'\$ reac \a, b) = (a, b\Q reac \a', b') are the amplitudes of reaction between reagents 
A and B (supposed to be known from quantum mechanics). We introduce the notation 
a,b and a', V for state of reagents A, B (atoms or molecules) and for state of atoms in the 
reaction products AB. Here, gj, q£, q£, q£ and q a ', Qw, Qa, Qb are the operators of creation 
and annihilation of atomic states a', b' for molecule AB, and a, b for A and B, respectively. 

Parameter of the reduced description are the averaged densities of gas atoms absorbed and 
not absorbed on the metal surface (n s (-R))' = Sp (^ha(R) Q(t)j , (h a (R)Y = Sp (^n a (R)g(t)^. 
ha{R) is the density operator for gas atoms adsorbed on the metal surface in state v\ 
ha{R) = ^2f a ^tj(fy^vj(R)i ^vjify are the creation and annihilation operators of 

gas atoms absorbed in state v on the metal surface which satisfy the Bose-type commutation 
relations. Since, we do not consider a catalyst surface explicitly in this model, the states v 
and /i mean the adsorption centers, where atoms can be located. h a {f) = Ylf" 5{f ' — r}) is 



the microscopic density of gas atoms, (G^(R, R')Y = Sp yO^(R, R')g(t)j is the nonequi- 
librium pair distribution function of atoms or molecules adsorbed on the metal surface, and 
Gll(R,R')=nZ(R)nl(R'). 

The relevant statistical operator 



9rel ~ Z R (t) 



-P \SH{t) dr^ a {r-t)8n a {r-t) (16) 



~ E E / "EE/ dRdR'M»£(R, R'; t)6%(R, R'; t) 

a v ab V V 

The partition function of the relevant statistical operator 



9-1 



Z R (t) = Sp 



q-l 



p(8H(t) " E/ dr^a(r-t)5h a (r-t) 



(17) 



"EE / / dm>M^(R,R>;t)G2(R,R>;t) 



ab 
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Sp(---) = Ua I W$0$mi s P{v,Z,<t){- ■ •)> = {-No, -Wo}, means averaged sum- 

mation over all values of spin and quantum numbers. The parameters /j a (r;t), /j^(R;t), 
M^(R, R'; t) should be determined from the corresponding self-consistency conditions 

(na^Y = (fiaWU K(R)Y = K(R))U (G U M R!)Y = (6%(R, R')Y rel , 

from where we fined that /i a (r; t) defines a local chemical potential of gas atoms; (J%(R; t) is 
a local chemical potential of an atom adsorbed in a state v on the metal surface. 5H(t) = 
H - (H)\ 5n a (f;t) = n a (f) - (h a {r))\ 5%(R;t) = h^R) - 5G^(R,R>;t) = 

G^(R,R')-(G^(R,R')Y. 

According to (3) and (17), NSO of the system "gas-adsorbate-metal" has the following 
form 

g(t) = g(t) rel + J2 J drf e £ ^T(M') (jf drgl el {t')I a {f- t) g 1 ^ (f) ^ a (f; f)) dt' 

EE/ di? f j^-^t') ^jf'dr^co^^OeJTCO^SK^o) dt' (is) 

EE/ dfttf?jT e^TW) (j\rgUtVG^ 



+ 



+ 



a,b v,n 



I a (f;t) = (l-P(t))^-\t)fi a (f), I»(R;t) = {l-P{t))-r\t)^{R), 
Q Q 

iZ^ti-A = (l-P(t))-1;-\t)G a (R,R!) 

ao q 



(19) 



are the generalized flows describing reaction-diffusion processes. The function ip{t) is defined 
by the relation 

q-l 



m = i 



-P I 5H(t) " E / d?ti a (r;t)6h a (r;t) 



(20) 



"EE/ *4«*)-EE / dRdR'M»£(R,R'-,t)G^(R : R'-,t) 

a v ab v)i , 

By means of NSO (19) we obtain the set of self-consistent generalized transport equations 
for averaged densities of adsorbed and not absorbed atoms along with nonequilibrium pair 
distribution function of absorbed atoms (or molecules). It has the following form 

+ EE I d& f e&-*^ftit;t,i!)prf(it;i!)tf 

+ EE / dn ' n " f e^V^^^ 



a'b Vp! 



^(K(R)Y = (^K(RlH'}Yrei + (^mR),H reac ])i el (22) 
+ EE 1*8 f ^"'^(a^O/W^O^ 

b v> J J -°° 

+ EE / d *'R" f ^V^gJ^^^)^^^)^, 

-,T i i J J —OO 



§- t (G^R, R')Y = (^[G^(R, R'), H'\)l el + (±.[G^(R, R'), H reac ])l el (23) 



y J J —oo 

+ EE / dR " f e £ ( ^ '-V^ ^ (i^,i^^/^'';^,^')^r( J R / ;^ / )^ 

-,T, ../../ J J —OO 



a'b' v'jx' 

where (f na n b i V 9 ™^' ^ V n a n- b i v'g^-G-,- are the generalized transport kernels. The second term 
in the right-hand side of (23), (^[n^(R), H reac }y rel , defines an averaged value of the operator 
of rate of the reaction between adsorbed atoms on the metal surface. Transport kernels 
are built on the generalized flows (20) taking into account contributions of amplitudes of 
chemical reactions in the flows I£(R; t') and Iq L _{R, R'; t'), and have the following form 

<Pbb> = Sp(l B (t)T(t,t') j\r e l el {t')I B ,{t')el7{t')y (24) 



In particular, (p nanb (r,r';t,t') describes dynamical correlations of diffusive flows of gas 
atoms and is connected to the inhomogeneous diffusion coefficient of atoms (or molecules) 
Dabi^r^'it). Similarly, the transport kernel ip™ n _(R, R';t,t') describes dynamical dissipa- 
tive correlations of diffusive flows of atoms adsorbed in states v and v' on the metal sur- 
face and determines inhomogeneous diffusion coefficient of atoms adsorbed on the metal 
surface D^' (R, R';t). Transport kernel ^_ nb {R 1 r';t,t') 1 (f, R'; t, t') describes dynam- 
ical dissipative correlations between flows of gas atoms and atoms adsorbed on the metal 
surface, and determines inhomogeneous coefficient of mutual diffusion "gas atom-adsorbed 
atom" D^(r,R';t). Investigation of these diffusion coefficients is very important. Trans- 
port kernel ^q_ iV {p — n , n) describes dissipative correlations of flows and densities of 
adsorbed atoms with flows of atoms, molecules and adsorbed atoms. Memory function 
(fr^lr. -(R, R' , R" , R'";t,t') describes react ion- diffusion processes between atoms adsorbed 
on the metal surface. They are higher memory functions with respect to dynamical vari- 
ables G-^. 

ab 

IV. CONCLUSIONS 

Consequently, we obtain generalized transport equations (22)-(24) for the nonequilibrium 
averaged densities of adsorbed and not adsorbed atoms of consistent description of the 
reaction-diffusion processes in the system "gas-adsorbate-metal" within Renyi statistics. At 
q = 1 these equations coincide with the equations of reaction-diffusion processes within Gibbs 



statistics 



851 ] . As we can see, these equations are nonlinear and spatially inhomogeneous. 



They can describe strong as well as weak nonequilibrium processes in the system. 
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